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We construct an approximate renormalization transforma- 
tion for Hamiltonian systems with three degrees of freedom 
in order to study the break-up of invariant tori with three 
incommensurate frequencies which belong to the cubic field 
Q(r), where + — 2r — 1 = 0. This renormalization has 
two fixed points : a stable one and a hyperbolic one with a 
codimension one stable manifold. We compute the associated 
critical exponents that characterize the universality class for 
the break-up of the invariant tori we consider. 



I. INTRODUCTION 

In this paper, we define an approximate renormaliza- 
tion scheme for Hamiltonians with three degrees of free- 
dom in order to study the break-up of invariant tori 
with the frequency vector = (t^ -t- r, r, 1) where 
T = 2 cos(27r/7) « 1.2469796 is the root of modulus larger 
than one of the polynomial : 

+ - 2r - 1 = 0. 

This approximate renormalization is defined for the fol- 
lowing family of Hamiltonians which are quadratic in the 
actions A = (^1,^12,^3) : 

H{A,v) = H,{A) + V{A,^), (1.1) 

where the Hamiltonian Hq is quadratic : 



Ho{A)=uJo-A + ^{n-Af, 



(1.2) 



and the perturbation is described by two scalar functions 
of the angles ip = {(pi,(p2,(p3) : 



V{A,^)=g{^)n-A + f{ip). 



(1.3) 



For Hq, the invariant torus with frequency vector is 
located at A such that fi ■ A = and ■ A = E where 
E is the total energy of the system. Since ujq satisfies a 
diophantine condition, the KAM theorem for Hamiltoni- 
ans (1.1) states that for a sufficiently small and smooth 
perturbation V, this invariant torus persists 0. For a 
sufficiently large perturbation, this invariant torus is bro- 
ken (converse KAM (|]). 

The purpose of this paper is to construct a renormaliza- 
tion transformation in order to investigate the properties 
of invariant tori with the frequency vector at critical- 
ity. The idea of the renormalization approach is to iter- 
ate a transformation in the space of Hamiltonians. For 
Hamiltonians that have a smooth invariant torus with 



frequency vector Wq, the iteration should converge to a 
trivial fixed point. The set of Hamiltonians that have 
a non-smooth invariant torus with this frequency vector 
form a surface (called critical surface) which is invariant 
under the renormalization, and is expected to be a codi- 
mension one stable manifold of a non-trivial fixed set of 
the renormalization. 

The aim is to define a renormalization as was done for 
Hamiltonians with two degrees of freedom in Refs. 
This renormalization is a combination of a partial elim- 
ination of the perturbation (we eliminate the Fourier 
modes of the perturbation which are sufficiently far from 
resonance) and a rescaling transformation which consists 
of a shift of the resonances, and a rescaling of time and 
of the actions. 

An approximate renormalization scheme for the fre- 
quency vector uio was defined in Ref. |^ following the 
renormalization defined for the spiral mean torus con- 
structed in Ref. Q . The main result of Ref. was that 
the renormalization has a fixed point on the critical sur- 
face but it is not hyperbolic. This renormalization dy- 
namics is structurally unstable. 

We propose to improve on Ref. just as was done for 
Ref. in Refs. |^,^ for the spiral mean torus. The main 
difference with the scheme constructed in Ref. ||^ is that 
we use a different normalization condition for the rescal- 
ing in the actions and we include a term that was previ- 
ously neglected. The main result we find is that a large 
part of the critical surface is the codimension-one stable 
manifold of a hyperbolic fixed point of this renormaliza- 
tion. We expect the dynamics of the exact renormaliza- 
tion to be the same as this approximate renormalization, 
at least locally. 

Following the approach of Escande and Doveil [|io|,0 for 
systems with two degrees of freedom, we perform two ap- 
proximations in this renormalization : 

(1) A three resonance approximation : we keep only the 
three main Fourier modes of the perturbation V at each 
step of the transformation. 

(2) A mean-value quadratic approximation : we neglect 
the dependence on the angles of the quadratic term in 
the actions. 

The frequency vector cjq is an eigenvector of the matrix 
N with the eigenvalue (r -|- 1)""'^ « 0.445, where N is the 
transposed matrix of the following matrix N with integer 
coefficients and determinant -1 : 



TV 




(1.4) 
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From N one generates a sequence of periodic orbits (with 
frequency vector {Jlfe}) approximating the motion with 
frequency vector ivq, i.e. (p(t) = uj^t + <^o niod 27r : 



J^fc = (Pk/rk^Qk/rk, 1) 



where pk, qk and rfc are determined by the foUowing re- 
cursion relations : 

Pk+i =Pk + '^qk + 
qk+i = Pk, 
Tk+i = qk + rk, 

and po = 3, go = 1, '^o = Ij i-e. ^k = N^'^CIq. These 
relations define a sequence of simultaneous rational ap- 
proximations {pk/rk,qk/rk) to the pair (r^ +t,t). 
The sequence of resonances is generated by the matrix 
from the vector i^i — (1, 0, 0) : 

The small denominators decrease geometrically to zero 
with the ratio (r -|- 1)"^ w 0.445 : 

u}o-Uk= T(r + 1)2-'=. 

The modes Uk are in resonance with the periodic motion 
with frequency vector Jlj, (with period r^) : 

^k ■ l^k+4 = 0. 

The matrix N can be decomposed in the following way : 
N^LPL^, (1.5) 

where 



L = 




and 



P 




(1.6) 



(1.7) 



We notice that det L — — det P — 1 and that P is a 
permutation of two elements of the basis; in particular, 
it satisfies = 1. This decomposition follows from the 
construction of the Farey sequence for an incommensu- 
rate vector a; G M'^ p^,[7|,p[ . With this decomposition 
of the matrix iV, we construct the renormalization from 
two operators (one associated with L and the other one 
to P) in a similar way as MacKay did for Hamiltonian 
systems with two degrees of freedom in Ref. |l^ ] in order 
to investigate the break-up of invariant tori with arbi- 
trary winding ratio. The approximate renormalization 
transformation we define for ujq can be constructed for a 
more general frequency vector from these two operators 
given the Farey decomposition of the frequency vector. 
For instance, concerning the spiral mean torus [HJ^, the 
renormalization is equal to one of the operators (the one 
which is associated with L) . 



II. DEFINITION OF THE RENORMALIZATION 
TRANSFORMATION 

The renormalization transformation TZ we define for 
the frequency vector u}q is composed of two operators : 
one associated with the matrix P and the other one with 
L. We denote TZp and TLl these operators. The transfor- 
mation TZ is defined for a fixed frequency vector ujq and 
is given by 

7^ = 7^L o7^L oTep oTeL. 



A. Definition of TZp 

The renormalization operator TZp acts on the following 
Hamiltonians : 



H{A,ip)^Ha{A) + h{n-A,ip), 



(2.1) 



where 



H^{A)=u,-A + -{n-Af, 



h{n ■ A,ip) = hi{A) cosip, = ^ifi + ■ A) cos 



It contains a shift of the Fourier modes and a rescaling of 
time and of the actions. The shift of the Fourier modes 
is constructed such that it exchanges the mode (1,0,0) 
and (0,0,1) without changing the mode (0,1,0). More 
precisely, we require that the new angles if' satisfy : 

cos ifi'i — cos ip3 , 
cos Lp'2 = cos ip2 , 
cos Lp'^ = COSipi. 

This is performed by the following linear canonical trans- 
formation : 

iA,v)^{A',v') = {PA,Pv), 

where we recall that P is symmetric and orthogonal. The 
vectors and Jl are changed into Pu: and Pi},. We im- 
pose that the images lj' and n' of the vectors u) and 
Jl by the renormalization TZp satisfy the following nor- 
malization conditions : fi' must be of (euclidean) norm 
one, and the third component of must be equal to 
one. Since P is orthogonal, PO is again of norm one, 
and thus O' = PJl. The third component of Pu) is 
equal to uJi. We rescale the time by a factor wi, i.e. 
we multiply the Hamiltonian by 1/wi. Then w' is given 
by uj' = u}/uji = (tJs/wi, W2/t^i, !)■ 

The quadratic part of the Hamiltonian Hq becomes 
{CI' ■ A)2/(2(jJi). In order that this quadratic term is 
equal to (O' • A)'^/2, we rescale the actions by a fac- 
tor Xp — l/i-^i, i-e. we change the Hamiltonian H into 
\pHiA/\p,<p). 
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In summary, a Hamiltonian H given by Eq. (2.1) is 
mapped into 

1 ^ 
H'{A, <f)=u:'-A+-{n'- Af + ^2(11 + 91^' ' A) cos ^„ 



where uj' = (wa/wi , 072/1.^1 , 1), ^2' = [fls, 512, ^1), and 





= hl^l, 




= h/^l 


fs 


= hl^l, 


9[ 


= .93/wl, 


1 

92 


= .92/^^1, 


1 

.93 


= .91/^1- 



The renormahzation operator 7?.p is equivalent to the 10- 
dimensional map given by the above equations (we recall 
that we impose a normalization condition on u) and Jl) : 

(u;, f,,g,; z = 1, 2, 3) ^ (a;', n\ gl; z = 1, 2, 3). 



B. Definition of TZl 

The renormahzation operator TZl is associate d to the 
matrix L and acts on the family of Hamiltonians It 
contains an elimination of the mode vi of the perturba- 
tion, and a rescaling procedure (shift of the resonances, 
rescaling of time and of the actions) such that the image 
of a Hamiltonian H given by Eq. (2J) is of the same 
general form as H and describes the system on a smaller 
scale in phase space and at a longer time scale. This 
operator is similar to the approximate renormahzation 
constructed for the spiral mean torus in Refs. 
We eliminate the mode Vi of the scalar functions / and g, 
by a near-identity canonical transformation. We perform 
a Lie transformation generated by a function 5(A, y). 
The image of a Hamiltonian H is given by : 



H' = exp{S)H ^H + {S,H} + ^{S, {S, H}} 



(2.2) 



where { , } denotes the Poisson bracket : 
{f,9} 



dip OA 



dg df 
dip dA' 



and the operator S acts on H like SH = {S,H}. The 
generating function is chosen linear in the actions and of 
the form : 



S{A, (p) = {z + yfl ■ A) sin (pi . 
We choose z and y in the following way : 

z = -/i/t^i, 

y = -{91 - fi^i/uji)/uji. 



(2.3) 



Then the generating function S satisfies : 

{S,Ho} = iQ{A)-hi{A))cosipi, 

with Q{A) = yfljjfl ■ A)'^. The image of a Hamiltonian 
H given by Eq. (2.1) is : 



H' = Hq + hi cos ipi + /i2 cos + ^3 cos ip3 

+ {S, Hq} + {S, hi C0S<y9i} + {5, /l2 C0S</J2} 

+{S, h3 COS ^3} + 1{S, {S, Ho}} + 0{e^). 

We neglect the term Q{A) cosipi produced by the terms 
{S,Hq} + hicosifi, i.e. we neglect the dependence on 
the angles of the quadratic terms in the actions. The 
term {S,hiCosipi} is of degree one in the actions and 
contains the modes and 2i/i. We neglect the mode 2i/i 
and we eliminate the mode of the linear term in the 
actions by a shift : A' = A + a where a is of order (in 
the direction of Jl). We neglect the modes with frequency 
vectors i^i ±1^2 produced by the term {S, /12 cos ^2}- The 
term {5, ^13 cos (^3} generates the modes 1^1 ± 1^3. We 
neglect the mode i>i + 1/3 and keep the next relevant 
Fourier mode i/i — 1/3 whose amplitude is denoted /13 = 
/3 -f .93$! • A, where 



/^ = (zg3f^i+2//3f^3)/2, 
g'^ = {Qi + ^(3)2/53/2. 



(2.4) 
(2.5) 



The term {S,Qcoaipi}/2 produced by {S,{S,Ho}}/2 
gives a contribution to the quadratic part of the Hamil- 
tonian. The new quadratic part is equal to m(0 • A)^/2 
where m is given by : 



m = 1 + -y'^nj. 



(2.6) 



Then after the elimination of the mode Ui, the new 
Hamiltonian is equal to : 

H' = uj ■ A + TO(r2 ■ A)'^/2 + /i2 cos ip2 + ^3 cos ips 

+h'3C0s{(pi - LPs), (2.7) 



where /13 is given by Eqs. ( |2.4D and (2.5). 

We shift the Fourier modes according to the following 

linear canonical transformation : 

{A,ip)^{A',ip') = {L-'A,Lip), 

such that cos (^2 = cosyj^, cos ips = cos ^^'^ cos((/3i — 
(fiz) — COS (^3. The vectors 17 and uj are changed in the 
following way : 

n' = = {^2,^3, ni - r!3)/(i + n^i^s - 201)), 
(^2/(^1-1), 1/(^1 -i),i). 



\Ln\ 

Lu) 



We rescale time in such a way that the linear term in the 
actions of H^ is equal to a;' : we multiply the Hamiltonian 
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H' by a factor 1 / (uoi — 1) . The quadratic term of the new 
Hamiltonian is equal to m||Zr2|p(r2' • A)^/(2wi - 2). In 
order to map this quadratic part into ($7' • A)^/2, we 
rescale the actions by a factor 

AL=TO||Lnf/(wi-l), 

by changing the Hamiltonian H' into XlH'{A/Xl,(p). 
After this elimination and rescaling procedures, a Hamil- 
tonian H is mapped into 

H" = a;'-A+i(ri'-A)^+/i"cos<^i + /i2 cos(^2 + ^3 cos(/?3, 
where h'^ = /f + g'/n' ■ A and 



/2 = /3C/, 
/•// f/ 

h = JsCf, 

9l = 52Cg, 
52 = 53Cg, 

93 = g'sCg, 

with c/ = m\\Lnf/{uji - 1)^ and Cg = - 1). 

The renormalization transformation is equivalent to the 
10-dimensional map 

(uj, 5,; z = 1, 2, 3) ^ (u;', n\ f^,^l- i = 1, 2, 3), 

given by the above equations. 

C. Renormalization transformation H. 

In summary, the renormalization TZ = TZj^TZpTZL acts 
in the following way : a canonical transformation elimi- 
nates the three main Fourier modes and produce the next 
three resonances which are the main Fourier modes at a 
smaller scale in phase space. A rescaling procedure shift 
these resonances into the original Fourier modes, and nor- 
malize the new Hamiltonian by rescaling the time and the 
actions. 

The rescaling procedure acts on in the following way : 
u}Q is changed into = (1/(3 - t^), l/(3r - t^), 1) by 

TZl, then into Wq^-* = (3 — T^,r~^, 1) by TZp, then into 
u}'^^ = (l/(2r - T^), 1/(2 - t2), 1) by TZl, and after the 

final step TZl, into o'g^'' — ujq. Thus the renormalization 
TZ is defined for a fixed frequency vector ijJq. In the same 
way, the action of 7^ on reduces to the map : 



n' = 



Nn 

\Nn\ 



(2.8) 



The spectrum of N is real and consists of the follow- 
ing eigenvalues -t « -1.2469, (t + 1)"^ «j 0.4450, 

the 



euclidean norm one) associated with the largest eigen- 
value of N. With this reduction, the renormalization 
TZ reduces to a 6-dimensional map {fi,gi]i = 1,2,3) 1— > 
ifi,gi;z = 1,2,3). 

The renormalization transformation we define can be also 
constructed for the following family of Hamiltonians : 

H{A, ip)=u;o-A+^A- MA + g{ip) ■ A + f{^), (2.9) 

where M is a 3 x 3 symmetric matrix with non-zero mean- 
value, and g is a three dimensional vector. Applying the 
renormalization changes the matrix M into : 



TZ{M) 



NMN 
tr (NMN) ' 



(2.10) 



1.8019. Thus by iterating the map (2 



vector fl converges to the eigenvector (denoted iV^, with 



where tr (NMN) is th e tr ace of the matrix NMN . 
By iterating the map ( ^.1C| ), the matrix converges to 
$7'^' ® ri^'^). Moreover, g is renormalized into Ng which 
tends to be a ligned to Jl'^' by iteration. Then the Hamil- 
tonians ( |2.9D tend to the degenerate Hamiltonians (1.1) 
under the iterations of renormalization. We expect that 
the Hamiltonians (2.9) belong to the same universality 
class as the Hamiltonians (|l.l[). 



III. RENORMALIZATION FLOW 

The renormalization transformation TZ has the follow- 
ing properties : TZ has an attractive integrable fixed point 
Hq given by 

i7o(A)=a;o-A+i(0(3).A)2, 

and another non-integrable fixed point iJ, which lies on 
the boundary of the domain of attraction of iJg. Outside 
the closure of the domain of attraction of Ho, the iter- 
ations of renormalization diverge to infinity. The renor- 
malization dynamics has the same qualitative features as 
the renormalization for the golden mean torus for Hamil- 
tonian systems with two degrees of freedom Q. The 
properties of critical tori are given by the analysis of the 
renormalization around the non-trivial fixed point ff^. 
In particular, the existence of a non-trivial fixed point 
for an exact renormalization would imply self-similarity 
of critical invariant tori. The stable manifold of iJ* is of 
codimension one; the linearized renormalization around 
-ff* has only one eigenvalue of modulus larger than one : 
5 « 3.4414. The value of the total rescaling coefficient in 
the actions is A « 11.2726, and the rescaling coefficient 
of time is r + 1 « 2.2469. 

We also find a critical fixed cycle with period 7 on the 
critical surface. This periodic cycle is obtained from the 
non-trivial fixed point iJ* by changing the sign of the 
Fourier coefficients /; and gi. The existence of this cycle 
is explained by symmetry reasons as it was done for the 
critical cycle with period 3 for the golden mean case H . 
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In particular, it involves the same critical exponents and 
scaling factors as the non-trivial fixed point, and belongs 
to the same universality class. 

IV. CONCLUSION 

We have defined two elementary operators TZp and TZl- 
With these operators, we have defined an approximate 
renormalization in order to study invariant tori with fre- 
quency vector ijJq = (r^ -I- r, r, 1). The renormalization 
has a hyperbolic fixed point with codimension one stable 
manifold. Consequently, we expect critical invariant tori 
with this frequency vector to be self-similar at criticality. 
In order to give a firm basis of this statement, it will be 
interesting to build an exact renormalization transforma- 
tion without the drastic approximations we used. 
We notice also that the approximate renormalization we 
define in this paper can be generalized to an arbitrary 
incommensurate frequency vector given by its Farey 
sequence, with the operators TZp and TZl- For a periodic 
Farey sequence the hyperbolic invariant sets are expected 
to be fixed points, periodic orbits or strange non-chaotic 
attractors (such as the spiral mean torus [[l5| ) . For a non- 
periodic Farey sequence, this fixed sets can be strange 
chaotic attractors, but this point has not been investi- 
gated yet. 
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